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LARGE DEFLECTION OF CURVED PLATES 
By H. G. Lew, J. A. Fox, and T. T.'Loo 

SUMMARY - 


Several problems on the large deflections of curved plates with com- 
pound curvature are treated. The large-deflection equations of plates and 
of cylinders are used. The effect of initial curvature for a plate under 
axial compression is to increase its deflection considerably upon the 
application of load. With increasing loads, the deflection curves merge 
into one regardless of the initial deflection form. It is shown also that 
the average shear strain for a plate under shear loading is not much 
affected by the initial curvature, at least for the initial deflection 
function used. In addition, for the circular cylindrical plate with small 
initial curvature under axial compression, the deflections of the plate 
are affected more by changes in the radius Rq than by changes in the 

initial curvature. 

It is shown that the effective width of the curved plates in longi- 
tudinal compression is reduced by the presence of an initial deflection 
function. 


INTRODUCTION 


The present report treats several problems on the large deflection 
of curved plates with compound curvature which are typical of aircraft 
structural elements. These elements are acted upon by edge loads (axial 
or shear) and the stresses and deflections are investigated for loads 
causing deflections of the order of the thickness of the plate elements. 

The work is divided into two parts. In the first part, the large 
deflection of a plate with small, initial curvature is investigated for 
two loading conditions: axial (compression) and shear. In the second 

part, the large deflection of a curved plate element with small initial 
curvature from a circular cylindrical shape is considered. 

In addition to the determination of the deflection and stresses the 
effective widths are calculated for several cases. 

A right-hand xyz coordinate system is used in this report. 
Rectangular-plan-form plates are considered and the origin of the coordi- 
nates is taken for each problem as given in figures 1, 2, and 3- 

This investigation was conducted under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
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SYMBOIS 


AqjBqjCrjDr coefficients given in appendix B 
a,b,l sheet dimensions 

Bi functions of w mn and 5^; used in equation (12) 

C coefficient used in equation (28) 


D = Eh?/l2(l 

d 

<ie 

E 

^■^qr^st 

F c 

Fo 

G 

h 


‘/ 8 


k = w o max /a 2 




unknown coefficients of complementary stress function (see 
eq. (13)) 

- V2) 

total, width of sheet, any coor dinat e system 
effective width of sheet (see fig. 7(a)) 
modulus of elasticity 

any stress function and stress-function coefficients 
complementary solution of equation (13) 
particular integral of F 
shear modulus 

thickness of shell or plate 
constants in equation (l6) 


M a . cosh ais + iWga 

1 2a sioh ass 

2a 
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^n,q,r,s,t 


max 

P 

Px 

^x-'Py 


*o 

U 

u,v/w 

u 

Uo 


constants 

maximum value; used as a subscript 
surface load 

load per unit length of edge 

average compressive stresses in x- and y-directions, 
respectively 

radius of circular cylindrical curved plate from which compound 
curved plate is a deviation 

strain energy of bending and stretching 

displacements in x-, y-, and z -directions, respectively 

average displacement 

initial displacement 


V total potential energy 

¥ potential energy of applied load 

VjVjjjn total deflection and coefficient of total deflection 

w n ,&mn initial deflection and coefficient of initial deflection 


v c deflection at center of plate 

v-j_ additional deflection 


x,y,z coordinates of a point 

a = w <W a 

6 ^ 1 } 6y 1 , 7xy ' median-fiber strains 

e x ",e y ",7 xy ." extreme-fiber bending and shearing strains 
average edge strain 


v 


Poisson's ratio 
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°xy 

cr * re 1 re 1 

°x >°y j°xy 

re U re U re " 

°x >°y y°xy 

n 

LARGE DEFLECTION OF CURVED PLATES 
Flat Plates With Snail Initial Curvature 


filer stress 
median-fiber stresses 

extreme-fiber bending and shearing stresses 
average value 


Equations of equilibrium .- The equations governing the large deflect 
tion of a thin curved plate in which the deflections are still small enough 
to use the simplified formulas for the curvatures have been derived by 
Von Karmen and by Marguerre (ref. l) and can be 'written as 

2 


V^F = E 


^2 
o v. 


>2 >2 
d a w^ 


-,2 

o w. 


^2 
o w^ 


[\dx V dx 2 V 


+ 2 


cix 8y cbc c)y 


2 2 


2 2 
d v 0 8 \ 


8x 8y Stx 


( 1 ) 


and 


v&v = e + a 

D D 1 


S 2 F 


( g i 






cbc* 


* 2( v ± + v o) 

$y 2 


- 2 


8 ^w-, 


+ v. 


8 2 f 

8x $y 8x dy 


ol 


( 2 ) 


Equation (l) represents the compatibility equations for the median-fiber 
strains for a plate with an initial middle surface deflection w Q . These 

strains are 


e y 


t 


'xy 


_ chi. , l/8w\ 2 _ if ^o ] 2 

8x 2\8 x/ 2\8x / 

_ 8v , i/dw\ 2 _ i /8w q \ ^ 

^y 2\cy/ 2\3y / 

_8u + 8v + 8w8w_ ^0 dw 0 
Sy cbc 8x 8y ck $y 


> ( 3 ) 
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Equation (2) is the equilibrium equation taMng into account the initial 
deflection w 0 of the middle surface. The additional deflection appears 
in the left side since it involves bending moments which depend only on 
the change in curvature while the slopes on the right depend on the total 
deflection. 


In addition, the median-fiber stresses and strains are 

0- * = SUL 
X 

a ' = 
y 8x 2 


, ' = - 
^ 8x 8y 

. _ l/b^F 


_ v 

EU2 . 21 

py ax 


. = l(~d 2 F _ dV 

V 8y 2, 


y ' = _P ^ + V ) ^ F 

/x y E bx. by 


The extreme-fiber bending and shearing stresses are 

a + v £"n 


c x = - 


- v 2 )' 


,8x £ 


°y 


a l d \ 


8 2 w. 


+ V 


2(1 - v 2 )W 2 bx 2 


.2 

Eh * W 1 


^ 2(l + v) by 




(4) 


(5) 
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The strain energy of stretching and tending of the plate is given hy 

U = ^ JJ j( ° x ' + < V') “ 2 ^ 1 + [fx'V “ ( T xy ) J | ^ + 


Note that 


8 (l) 2 w 1 
dx 2 8y 2 


3 %. Y 


dx dy = 0 


if w-|_ is zero at the edges. Thus, the total potential energy is 

V = U + W (T) 

where W is the potential energy of the applied loads. The multiple 
integration is taken over the area of the plate. 

Plates with normal edge loads .- The coordinate system used for and 
the manne r of loading plates with normal edge loads are shown in figure 1 . 

Since this theory will he applied to plates with symmetric ini tial 
deflections and loads the assumption that there is symmetry about the 
center lines of the plate will he used here. The following boundary 
conditions appear to he reasonable representations of such plate elements: 

(l) The resultant load is constant in the x-direction (zero load 

in y-direction) ; that is. 



dy = -p^hb = Constant 


where P x is the load per unit length of the edge. 

(2) The edges x = 0 and x = a are straight. 

( 3 ) Shearing stresses are zero at all edges. 

For this case the potential energy of the applied load for normal 
edge loads along the edges x = 0 and x = a is 
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W 


= -/(°*' u )x=0,a ^ 


( 8 ) 


It is convenient, “with the coordinate system shown in figure 1, 
to assume v and w Q as a double sine series in x and y. Hence, let 


w = 


%m sin sin ^ 


m=l,2 } 3 ) — n=l,2,3,... 


(9) 


w D = 51 Bnm. sin “ sin 222 (lO) 

n=l,2,3.>»*» 

Upon substitution of equations (9) and (lO) into equation ( l) , it is 
seen that a particular integral for F will he 

2: 2: v c 0S as cos js (11) 

d q=0,l,2,... r=0,l,2,... a b 

where p^ is the average compressive stress in the x-direction. The 
value of Fqp in terms of the unknown coefficients w mn and S^m is 


00 



(12) 


Values of F^ r are tabulated in table 1 and the quantities are 

functions of w^ and 8^. 


The requirement that only one pair of edges x = 0 and x = a is 
loaded and that the other pair is stress free adds a complementary solu- 
tion to F (eq. (ll)). This function is 
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Assertion of the boundary condition 

°y - 0 

at y = 0,b leads to the value for of 


(HO 


c a = 


q.=2,lj-,6, . . . 


■qr 


cosh 


qrtb 

qrtb 2a 


2a 


sinh 


qrib 

2a 


(15) 


Thus, all the boundary conditions are satisfied by these solutions. 
Condition (l) is satisfied by using equations (3).> (ll), and-(H4-). Con- 
dition (2) can be verified if the relation between strain and displace- 
ment u is integrated. Th is integration leads to 


u x=a - **' 


' u x=0 


f a . 

£ - _ 1 

■v 

/W\ 2 _ (^o\ 2 

1 a/2 

x 2 

\8x/ \8x / 

4 


?dx 


= _ JLfe) - *!_ £ Z 

^ ; l£a . . . n=l,375, . . . V ” m 



which is a constant. The remaining condition of zero shearing stresses 
is obviously satisfied. 

Method of Solution .- The direct method of solution utilizing the 
principle of minimum potential energy will be -used instead of attempting 
to satisfy the remaining differential equation. Thus far, the stress 
function F is known in terms of unknown coefficients of the additional 
deflection w-j_(x, y) . By the direct method, a first variation of the 

total potential energy set equal to zero will yield nonlinear algebraic 
equations for the coefficient Wmn which then may be solved conveniently 

by successive approximation. 

The formation of the total potential from equation (7) by insertion 
of equations (9), (lO), (ll), and (13) leads to integrals which are listed 
in appendix A. 

In this procedure a system of nonlinear algebraic equations is 
obtained for the coefficients w^ in terms of the parameter of the 
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load (hp^h^/ jt%)) . For the numerical examples given, four unknown coef- 
ficients v^i, w.^., w^, and are used. It is convenient in the 

solution of these algebraic equations to assume a value for w-, and 

then to transpose all nonlinear terms to the right side. A first approx- 
imation is then obtained by neglecting all nonlinear terms on the right 
side and solving the remaining equations by Crout' s method. A second 
approximation is obtained by substituting the results of the first 
approximation in the right side of the equations and applyin g Grout's 
method. The convergence is very rapid. An example of this is shown in 
table 2. 


The numerical results indicate that three to four coefficients w mn 
are sufficient to represent the complete solution for practical applica- 
tions. In table 2, the center deflections w c /h using the four coeffi- 
cients v in , w n w 3i> anj ^ 8X6 gi' veri * There is also shown in 


*11' w 13> 


this table the coefficient 


33 

v n.l h 


which is equal to w c /h if only the 


w n i coefficient is used (with w-^, w^-p and w^ taken equal to zero). 
Thus it is seen from table 2 that the dominant coefficient is w-| j , and 


in most cases the fourth coefficient w^ is a very small percentage of 

the first. This is expected since the final deflection is assumed to be 
of the same form as the initial deflection. That is, the plate will bend 
more natur ally to the given initial deflection by a compressive load for 
small deflections. 


Numerical examples are considered for two maximum' initial deflec- 
tion coefficients of h and (l/3)h where h is the thickness of the 
plate. Comparison with the result from reference 2 is good if account 
is taken that a maximum initial deflection coefficient of O.lh was 
used there. The deflection curves with increasing load are shown in 
figure 4 together with the results from linear theory. This figure 
shows the coincidence of results (small and large deflections) in the 
linear range. Moreover, all deflection curves merge for large deflec- 
tions regardless of the different initial deflection amplitudes. Of 
course, this is physically evident; also, the merging point appears to 
occur outside the range of validity of the equations used. 

This method of solution (Ritz method) can be easily applied to 
large-deflection problems of a similar nature. Further calculation, 
for example, can be performed for initially curved plates with the same 
boundary conditions as those treated above from the integrals given in 
appendix A. In summary, it should be mentioned again that one differ- 
ential equation (for F) and the boundary conditions are satisfied exactly 
and the remaining differential equation is satisfied on the average by 
replacing it by the assertion of the minimum of the potential energy. 
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Plates under edge shear loads .- In this example, a square plate 
■with a small initial warp is considered. The plate is loaded in pure 
shear along the edges which are simply supported. The coordinate system 
is shown in figure 2. The warping is again treated as the initial 
deflection w Q of a flat plate. Thus, with a linear variation along 
the lines of constant x or y, the initial deflection is 

w Q = k-jjx + kgjr + kxy (l6) 

It Is apparent from equation (l) that the linear x and y terms of 
the initial deflection equation (l6) cannot have any influence on the 
stress function F: furthermore, these linear terms do not enter strain- 
energy expression (6) which is the same for this example. Thus, the 
coefficient k of the xy term in expression (l6) is given hy 

w 0 = kxy = ■°^ c xy = | xy (17) 

a 2 a 

where a = w<-, mBx j s '* ' w 'o mBX is as indicated in figure 2. 

In addition to simply supported edges, the plate elements will 
satisfy the following houndary conditions which appear reasonable for 
physical applications. At x = 0,a, ff x ’ = 0 and 1 = a^, and, 

at y = 0,a, Oy 1 = 0 and o^T = o^. It is noted that these condi- 
tions are satisfied, on the average. 

The deflections w^ under load are assumed as 


v-> = 


mn=l,2,3, • • • 


w™ sin 225= sin HZS l 

rrm a a 


(18) 


From symmetry considerations, v mn = % for all values of n and m, 
and, if m + n is an odd number, the terms must vanish for the same 
reason. In the example chosen, the series is truncated at n = m = 3 
since the leading terms dominate. The deflection w = W]_ + w Q is then 

given by the expression 

£ = sin — sin 2SC. + ^22(sin 2S sin 22SL + s in sin + 

£l|_Ji a d n \ a a a a / 

— sin 2ioc sln + ^22. sin 2l2= S in + . . .1 + 
h a b h a b J 

<m) 


(19) 
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The substitution of this expression into equation (l) yields a solution 
for the stress function F in the form 

Fgj. cos 2 212E. cob 2 2^ (20 ) 

* ^0 r=0 a 


The values of F are given in table 3- 

If equations (19) and (20) are substituted into strain energy 
expression ( 6 ), the strain energy U may be expressed in terms of the 
unknown parameters w-^, w-^, M 22’ an “ i w 33* ® ie potential energy W 

of the applied shearing forces along the edges is 



where the integration is performed over the plan form of the plate. The 

expression for ( + ^ ) can be found from equation (3) * 

\qy 

The total potential energy V (eq, (7) ) is now made stationary with 
respect to each of the unknown parameters. The solution of the four non- 
linear algebraic equations follows the iterative procedure of reference 3- 
The results are given in table 4 and figure 5- 


The problem has been evaluated numerically for V n = 0 and 
Wpmax = ® ie relationship between ^ ~G^")(h)^ + ^)(h)^ 


in each case has been plotted in figure 5* It is noted that, in the 
first case of an initially flat plate, the postbuckling strength of the 
plate increases continuously after the critical buc kling load is reached, 
whereas the classical small-deflection theory predicts a constant strength 
(broken line in fig. 5)* 


The stress field developed in the plate as the buckling proceeds may 
be readily computed from the values of the parameters w-^, w-^, w pp , 

and w^2 "by using equations (4) and (5) • 

In figure the results of an approximate solution (ref. 3) of 
Von Karman's large -deflection plate equations are plotted. This solution 
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assumes the form of W]_ identical with equation (19), except for a con- 
stant factor. However, scmewhat more rigid edge conditions were imposed 
on the plate (i.e., e x ’ =0 at y = 0,a and =0 at x = 0,a). 

This accounts for the slightly larger value of o^. for a given value 
of average shear strain. It is noted that this difference is small so 
that it appears that the edge conditions have a small effect on the rela- 
tion between applied average shear stress and average shear strain. 


Nearly Cylindrically Curved Plates 

The small-deflection equations of a curved plate formed from a body 
of revolution under edge and surface loads have been considered from the 
equations of equilibrium of a body of revolution ih reference 4. In this 
section, the large deflection of curved plates with small compound curva- 
ture is investigated by modifying the large-deflection equations of a 
circular cylinder to include a small initial curvature. The coordinate 
system and an Indication of the plate shape are given in figure 3 • 


The following equations govern the deflection of a curved plate 
which differs slightly from a circular cylindrical shape. It is noted 
that these equations immediately reduce to those of reference 5 (where 
other references are given) for the case of w 0 equal to zero. 


The median-fiber strains are 



_ 8u + 

) 2 + ^i 

^0 


/ 

dx 

Sbc 



2 



= & + 1&. 



y 

$y 2 \dy J 

1 Ro 


y 1 

_ chi , . 

dWf dw^ 

8w 0 

7 xy 

$y dx 

dx $y 

dx 


The median-fiber stresses are given by 

2 


a ' = E_ 

°x 


8u l(^l\ , ^1 


rt7 S + S\S7 


+ dx ax + v 


tr + 2 W 


Wi 

Ro 


« = _J 

y 1 - v 2 


J^ + l&) _^ + 

1 $ 7 2 W 1 Ro 


. Af + 

&x ^ \cbc / 8x 8x 


(22) 


J xy 


I __ 


2(l + v) \$y cbc . ck cfy- dx 8y / 


E 


Y (23) 
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In addition to the above stress and strain expressions the equili- 
brium equations in the x- and y-directions are, as before. 


8cr„ 


dcr. 


8x 


JL-0 

py 


Sow' Scfy' 


( 24 ) 

(25) 


If one introduces a stress function F defined by the set of 
expressions 


°x = 




CT » = 


°xy ' 


8 2 F 
8 x cty 


(26) 


and combines these with equations (22) through ( 24 ) a compatibility 
equation determined as 


vV = 


E 


( 

8 2 w 

-\2 
0 w 

1 8 2 w 

2 

3 w o 

\8x 87/ 

8 x 2 

$y 2 " 

R o a* 2 

ctx 2 8y 2 


(27) 


results . 

In place of a third equilibrium equation relating to the sum of the 
forces in the normal direction to the surface, the potential energy of 
the system is used for the direct method of solution as in the previous 
cases. One notes that a first variation of the total potential will lead 
to the equilibrium equation normal to the surface. The median-fiber 
stresses and strains and the extreme— fiber bending and shearing stresses 
are the same as those given by equations (4) and (5) if it is remembered 
that the coordinates x and y here are the cylindrical-surface 
coordinates . 

The strain-energy expression involving the sheet bending and 
stretching is given by equation (6). However, the boundary co ndi tions 
imposed on this problem necessitate another stress function. The poten- 
tial energy of the edge load will be the same as in the previous section. 
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namely, the average edge load, times the displacement. Since the edge 
displacement is a constant, this amounts to expression (8) with the 
boundaries x = ±1 substituted for x = 0,a. The initial deflection 
from the circular cylindrical shape will be taken as 

w 0 = Cx 2 (28) 

indicating a deviation of the meridian curve from straightness. 

Clamped conditions with regard, to the slopes at the bulkhead, or 
stringer junctures are assumed. The reason for this assumption was that 
it was felt that many times the bulkheads or ribs are fairly stiff in 
torsion about their flanges and that the flanges are quite wide, thus 
presenting a nearly clamped condition to the sheet ends. Furthermore, 
since the plate is initially bulged along its axis, it would probably 
continue to bulge under load which would mean that the sheet would act 
as a clamped sheet across stringers as well. Obviously, other cases 
could be considered where nearly simple support conditions might exist. 
Frcm the nature of the problem just described, it is necessary that the 
edges associated with the stringers move in radial lines j that is, the 
displacement v at the edges y = ±b must be zero. Also, it is assumed 
that the edges represented by bulkheads or ribs remain straight; that is, 
the displacements u at x = ±1 are constant. 

To satisfy the clamped condition the additional deflection function 
was assumed as thd series 

, Vl - i + °o» ®?) (l + oos Sf ) (29) 

Substitution of equation (29) into equation (27) will yield an expression 
for a particular integral for F. It will be 



where p„ and .p are the magnitudes of the average compressive stresses 

x y 

in the x- and y-directions, respectively. The values of F g ^- are given 
in table 5 and are similar in makeup to the other problems. 

Assertion of these boundary conditions adds a complementary solution 
to F which has- the form 
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F = E cob cosh Bq sinh + 

q=l i \ Z Z Z / 

E J ^(Cr c°si StSl) r «Uia) (31) 

The values of Ag, Bq, C*-, and Dp are given in appendix B. After 
applying boundary conditions to the complementary solution for the purpose 
of determining the unkn own coefficients, the displacement u is 


U ±l 


+ 


x 

4 



m=l n=l r=l 


Vmr m 2 * 2 - 





(32) 

The expressions for the nonzero integrals of the total strain energy 
in terms of the undetermined deflection coefficients are given in appen- 
dix C. Again, a solution vas obtained by collecting linear terms and 
first omitting all nonlinear terms except one which was specified. The 
equations were solved by Crout's method and the solutions obtained were 
used to evaluate the nonlinear terms which were omitted. This iterative 
process was carried out three or four times to yield three-place accuracy. 
It was most advantageous to specify w^ rather than w-q_ in this exam- 
ple for rapid convergence. 

Three numerical examples were chosen for this section. One plate 
has radius Rq equal to one-half the value of the other two; another 
has a maximum initial deflection of one-half the value of the other two 
and the same radius Rq as the first case. The results indicate that 
the effect of initial curvature is very small on the additional deflec- 
tion as compared with the effect of changing the radius Rq. The maximum 
additional deflections are tabulated in table 6 and the stresses are 
given in table 7 . The maximum additional deflections for one, two, and 
four terms are given in table 8 . 

A plot of the additional deflection at the center is shown in fig- 
ure 6 for the three cases considered and demonstrates graphically the 
effect of cylinder radius and initial curvature. 

The ratio of effective width to initial width defined as the ratio 
of the load carried by the plate to the load the plate would have carried 
if the stress had been uniform and equal to Young's modulus times the 
average edge strain (ref. 6 ) is shown in figure 7 . An equivalent 
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definition of the effective width in many cases is the ratio of the 
average stress across the loaded edges to the maximum stress at the edges. 
The results obtained using this definition are also shown in figure 7( a ) 
and clearly indicate the difference of the two definitions in the cases 
considered here. The reason for this difference is the assumption of 
straight edges with zero displacement of the y edges generating an 
average load p^. 


The effect of an initial deflection function for the curved plates 
loaded in longitudinal compression is a reduction in the effective width. 
In fact, the sheet deflects immediately upon loading since there is no 
critical buckling load. As the load increases, the outward (negative 
Z) bulging of the nearly straight circular cylinder Indicates a higher 
value of effective width than that obtained with the clamped flat plate 
of reference 7} this is probably due to the large curvature of the basic 
cylinder as compared with that of a flat plate. This is shown in fig- 
ure 7(a) . Figure 7(^>) compares the effective width of a single finite- 
length wave deflection of the problem considered here with the simply 
supported very long circular cylinder of reference 6. The increase in 
effective width here is due to the clamped boundary and the outward 
bulging assumed. It is noted from the figures also that the change of 
effective width with increasing edge strain is small. 


CONCLUDING REMARKS 


The problem of the large deflection of plates with compound curvature 
has been considered in this report. The large-deflection plate theory 
and circular -cylinder theory with an initial nonzero deflection function 
are used. Three cases have been worked out completely. These are an 
initially curved plate subjected to a longitudinal compressive load and 
to a shear load and a curved plate, with an initial deflection from a 
circular cylindrical shape, loaded with longitudinal compressive loads. 

The solution by the Ritz method i n dicates that a very small number 
of unknown coefficients (two to four, say) was sufficient to give reason- 
able results even with large deflections. Thus, the stress conditions 
and large deflections of curved plates (compound curvature) may be easily 
obtained by such methods. Same results are shown in the figures and 
tables. 

The effect of an initial deflection function for the curved plates 
loaded in longitudinal compression is a reduction in the effective width. 
The effective width of the curved plate with almost circular cylindrical 
shape is higher than that of the initially flat plate and that given by 
Levy (ref. 6) for circular cylindrical curved panels. 


The P enns ylv ani a State University, 

University Park, Pa., April 29, 1 955* 
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APPENDIX A 

INTEGRALS OP RECTANGULAR PLATE LOADED IN EDGE COMPRESSION 


The deflections of a rectangular plate loaded in edge compression 
are given "by the following expressions: 




m=l n=l 




ruty 


m=l n=l 


sin . 
a h 


If Poisson's ratio is taken as 0.3 for the material together with the 
deflections given above, the nonzero integrals are: 


r K' + =v') 2 tedy=i f~ t n sLjVi. !. 

esJ 0 J 0 {* °y) A 2 / ar 


h_ 

2E 


,3„3 sinh 


qjfb 


" 2a 


q.=2,4,6 r=2,4,6 t=0,2,4,6 a 2 


^qr^qt 


, . 03 4 4 

‘- 1 v rit ^p^. 2 + ^ «:sb F 2i 


" 

22 V*4j?,6 ■?“ “ '<* * q=lj,6 T T V 


■f- qV sinh 

4^,6 ^,6 a 2 «q VqO 


2h 

E 


_3„3 ^ + sinh S2* 
Jr a a 


d=2,4,6 r =0,2,4, 6 t=0,2,4,6 4 a 2 


M 2 . F l rP Q-i + 2E P x ab 



QlCVI 
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h£ r * *>kv - * * - 


0 u 0 


00 00 


^i(2.6 ^ ^ ^ — — 2“ a 

q=2,4.,6 r=0,2,4,6 t=0,2,4,6 4a 


p p p 

1 p t> + -u 2 qrtb 

" h 1 - —_P' C °^ 5T 


v 


■^qr^qt 


r r & a(w - m ] 2 «* * - ^ m g 5 n j 5 j p + “ 2 $) ow - v) : 


00 00 


dt " " §P^> - * 2 £ m2 ( w mn 2 “ 5 mn 2 ) 


where 


Mq. = cosh + 


qrtb 

®* . 2a 


sinh 


qrtb 

2a 
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APPENDIX B 

EXPRESSIONS FOR Aq> Bg, C r , AND D*. FOR CURVED SHEET 
HAVING 1 /Rq MUCH GREATER THAN MERIDIAN CURVATURE 


When 1/Rq is much greater than the meridian curvature and with 
Ujjy* = 0 along all of the edges , 

\ - M 1 + ¥ «*» 4 ) 

Cp = -D^l + ooth SSij 

Substituting these values of and C*. and assuming constant u 

displacement and zero v displacement will give the following values 
of and D^: 


<1 


1 

k 1 9 toht57h ( V U 2 + V U V 51 + + ^13 V 3j) - £h®o 

2 

S imE-f5E7Ty( v Il 2 + ^13 2 + a "U*31 + ^ 13 ^ 3 ) 

5 

w 1 sshtey( v u- v 3i + ^33 + V + ^33 2 ) - ^ LS iM 

k 

£ «lnh ?to/l)( W U V 31 + 9V 15 V }3) 

6 

T3 35 h BinhT^b/l)( V 51 + *"33 ) 

r 


I 

& l orfT7S)(^ 2 + ^ + ^31 2 + *W») - ™ 

2 

512 ! ,-5h~(ari/b)( v a a + •WiB + + 

3 

il i J,, + v 2 + 0v v t ^ 2 \ CU> ( V 13 + “33 ) 

192 1 Binh fetl/bA 11 13 V 13 ^5133 ^33 y olnh (J*l/b) 

k 

512 i olnh ?4s2/b)( V li V 13 + 9 u 31 v 3;) 

6 

1^5 7 Blnh (Ssl/b)( v 13 + ^33 ) 
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APPENDIX C 


INTECSAIS OF RECTANGULAR DOUBET CURVED PLATE WITH 1/Rq 
MUCH LESS THAN MERIDIAN CURVATURE 


The deflections of the rectangular doubly curved plate with 1 /Rq 
much less than meridian curvature are defined by the following expressions 


W 1 = 


2 


m=l,3.,5.» • • • • • - 


+ cos 5fflE)(l + cos 


w Q = Cx 2 


With these deflections^ the nonzero integrals are 

*££ * v)2 ^ ^ '•**(* * * s2 *) 2 + 

^ F * ■“ ¥ 


Whi ^ 

r=l,2,3.> • • • 


q 2 r 2 

^ iV 


^(-i) r sum ^ 
rrt l 


2V?v^v, 


£(_1)1 slnh «» 
C£TC b 


-A £ 

- - - t* 



PICVJ 
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' feXiXb 2(1 + v) [>V " (V) 2 ] ay = - Y (1 + 


where 


\jr “ 


qrZb 


2 , 2 . 2,2 

q. d + r Z 


£ Jl AA - * *A$ + + v Af + A 


V 2 (2- + i\ + * 2 /’9- 
*■ u 2 h 2 ; 55 lz 2 


+ 4Y 


h 2 ; 


+ 


k 

JC b f 2 


A 2 * ^ll W 15 + 81 v 31 2 + w 13 2 + ^3^35 + 8llf 33 a ) 

4 / v 

*3 VU. 2 + a, ll w 5i + 81 w 13 2 + i 6 sw 13 w 33 + w 3i 2 + 8lw 33 / 


8 Z 

4 

8 b 



22 


NACA TN 36^4 


REF EREU UES 


1. Mar guerre, K. : Zur Theorie der gekrupnten Platte grosser Formanderung . 

Proc. Fifth Int. Cong. Appl. Meeh. (Sept. 1938, Cambridge, Mass), 
John Wiley &■ Sons, Inc., 193 9 , pp. 93-101. 

2. Coan, J. M.: Large-Deflection Theory for Plates With Small Initial 

Curvature Loaded in Edge Compression. Jour. Appl. Medh., vol. 18, 
no. 2, June 1951* PP* 143-151 • 

3* Levy, Samuel, Fienup, Kenneth L., an d Woolley, Ruth M.: Analysis of 

Square Shear Weh Above Buckling Load. NACA TN 9^2, 1945* 

4. Lev, H. G. : Bending of Thin Plates With Compound Curvature. NACA 

TN 2782, 1952. 

5. Kempner, Joseph: Postbuckling Behavior of Axially Compressed Circular 

Cylindrical Shells. Jour. Aero. Sci., vol. 21, no. 5* May 1954, 
pp. 329-335. 

6. Levy, Samuel: Large-Deflection T h eory of Curved Sheet. NACA TN 895* 

1943. 

7. Levy, Samuel, and Krupen, Philip: Large-Deflection Theory for End 

Compression of Long Rectangular Plates Rigidly Clamped Along Two 
Edges. NACA TN 884, 1943- 



NACA !SN 3684 


TABLE 1 


VALUES OF F qr FOR w = 


• • • 

v Q = Bn sin ^ gin 2S[ AED a /b = 1 


Z % Bin S sin SSE 


(E/52)( W Z - e n 2 " ^11^15 + 9W31 2 " ^31X33) 


w n - 

1 

co 

w n v i3 

+ 9w 31 i 

'( V W 2 + 

v ll 2 - 

e 2 

5 11 “ 

w n w 3i 

+ v ll w : 


(e/400)^-v i;l v 13 + 9w i;l w 55 + 25 w 15 w 5l ) 

(e/400) ( 9 ^ 23 ^ 33 ) 

(e/64)^w 11 w 31 + 9w 13 v 33 ) 
(e/1)-00)^-w 11 w 31 + 9w 11 w 33 + 25 v 13 v 31 ) 
(E M)(-v 15 v 31 ) 

(e/ 2, 70^ ) (-9^13^33) 

(E/ 288 ) ^w 31 2 + 9w 33 2 ) 

(eAoo)(9w 31 v 33 ) 

(e/2,704)(-9w 31 w 33 ) 
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YAIUES CF F 

TABLE! 5 

FOR 8 SPARE, THTWAT.Tr CURVED PLATE UNDER SHEAR LOADING 


(l/ 8 )^!! 2 - 2»n»i 3 - ^13*33 + 9WJ3 2 ) 
+ 9 w 15 v 35 + 2wg2 2 ) 

(1/72) (v 13 2 + 9 w 35 2 ) 

( 1 A ) ( V H W X3 “ + (l/ 1 ^) w 22 K 

(l/200)(9v 11 w 55 - VjjVjj + 25^ 15 2 ) 

(1/200) (9^v 53 ) 

( 1/200 )(9v 11 v 33 - + 25v 13 2 ) 

(1/200)^5553) 

(1/1,552) (-9 v 13 v 33 ) 

(1/1,552) (-^^33) 

(1/52) (-v 13 2 ) 

^H" V 13 V 22) + 0A>«U* 

(l/50)^*VjjVgg + l6v~i 3 'tf 90 ) + (5/l00 

(l/55®) ^^13^22 "** ^ ^22^55) 

(l/^otfto^Vgg + 16^x3^22) + ( 5 /lOO)Kw 13 
(1/558)^16^^22 + 56 w 22 V 33 ) 

(1/578) ^WjjWga) 

(1/578) (-^ V 13 V 22) 

(l/56)w 33 K 
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TABLE 4 

NUMERICAL RESULTS FOR INITIALLY' CURVED SQUARE PLATE UNDER SHEAR 
LOADING FOR DIFFERENT INITIAL DEFLECTIONS 


(a) w D = 0 


V 22/h 

/ , .. \ 

V 8 - 2 

P 

V lc/h 

[ du + 6v\/a\ 2 

°xy a '" 

ISy cfcc/v 1 / 

2 

Eh 

Gh 2 

0 

22.087 

8.495 

22.087 

0 

.05 

22.857 

8.747 

22.745 

.592 

.ID 

24.612 

9.550 

24.257 

.759 

.15 

26.985 

10.091 

26.236 

1.037 

.20 

29.698 

10.959 

28.441 

1.295 

• 25 

52.659 

11.842 

50.789 

1.526 

.50 

35.854 

12.795 

35.268 

1.757 

.40 

42.921 

14.839 

38.582 

2.123 

.50 

51.005 

17.IOO 

44.461 

2.479 

.70 

70.692 

22.405 

58.255 

5.157 

1.00 

110.470 

52.770 

85.202 

4.159 


Cb) v D = 8 h 


CVJ 

CVI 

> 


O 

p 

v lc/h 

( du + ^V a f 

V 


V^ry- l&)\h) 

_2 

Eh 

2 

Gh 

0.05 

14.349 

4.401 

11.443 

1.003 

.10 

21.897 

6.647 

17.282 

1.438 

.15 

28.324 

8.505 

22.uk 

1.777 

.50 

47.583 

13.966 

36.311 

2.615 

.50 

70.873 

20.168 

52.438 

5.595 












28 


MCA TW 368^ 



[vts* ♦#](-& 

[yi?*i)!][-£-: 


v ll v 33 “ 16 v 13 v 31 


f v 15 v 53) 


-sw 


r ll v 31 " x6 V H V 33 ■ 16 V 31 ■ s W 3l" 


32 v H v 3l - 32 v 3l "52 V 13 W 31 “ 32 w 3l w 33 " % 


*33*31 "l *31*33 ‘55 "11*33) 


v u w 31 ' 32 v u v 33 - ’WSl) 


[vt? *■ ?)] ■ § *13 v 31 • S v 13”33 - § v 3l’33 ~ tZ "33 

[/(? * (' 6 "»*» ' 3§ WS - X "11*33) 

[vt 1 ? + s)3 (- § V U- ¥ 31 - 5 *U*33 - 3 a v «"3i) 

[“/Pf* + ¥)](■ b *u*53 ■ S *13*31 ■ 4 *u*3i) 

a al) 5 ] (- g - n - 5J - g - g w) 

[yt^*^)](-W33) 

[vf 3 ! 6 * s)3(- b V ‘ b "31*33) 

Vt 2 ? *¥)!|(-! *31*33) 
j/(^ 8-31*33 -S*33 2 ) 
G/( 3 f* 1 ! 1 )](-S *31*33) 


'*(£*5$ 


-2-v 2 _ JL v v ^ 

32 31 32 31 V 33J 

' £ v 3i v 33) 
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TABLE 6 

DEFLECTIONS FOR NORMAL LOADING ON CURVED PLATE WHERE Z/b = 2 
AND 1 /Rq IS MUCH GREATER THAN MERIDIAN CURVATURE 


w 3 V h 


-0.25 

-.50 

-•75 

- 1.00 

-1.25 

-1.50 


(a) Rq = 0.4b 2 /jth 


w 0 = h/ 2 


P x 1 2 /eA 2 v lc ^h p x b 2 /at 2 h 2 w lc/ /b 


O.689 

1.244 

1.734 

2.252 

2.858 

3-584 


-0.138 

-.251 

-.346 

-.438 

-.531 

-.636 


1.247 

-0.247 

2.242 

-.431 

3.566 

-.616 

1 


(b) R 0 = 0 . 2b 2 / jth j w 0 


p x b^/lbt 2 b 2 


0.553 

1.061 

1.529 

1.971 

2-392 

2.802 


w lc/ h 


-O.O58 

-.112 

-.163 

-.210 

-.256 

-.300 
























TABLE 7 

HCMERICAL VALUES OT BTBBB8BB CALCULATE) FOR DOUBEf CURVED SHEET WCTH if b a 2, 
1 /Rq MUCH ORBATER THAU HBEUDIAH OUKVAIUKE, AHD h 0 » h 




x position 



0 

tl 

0 

0 

tl 

0 

ii 


p? 2 
*2 2 

y position 

Ro ■ 

0 

0 

t b 

0 

0 

ib 

ib 


EX n 

a "b 2 

X 

eA 2 

bA 2 

°A 2 

■A® 

rf 

5 lL 

bA 2 

E A 2 

Y^- 

bA 2 

CVJ OJ 

of 
b « g 

OL't 2 

bA 2 

a ’b 2 
x 

bA 2 

,2 

eA 2 

a ’b 2 

X 

lA 2 

EjtV 

®x'b 2 

bA 2 

45! 

eA 2 


0.353 

• 0*092 

-0.088 

O.O72 

0.023 

0.020 

0.065 

-0.549 

0.098 

-0.542 

-0.250 

- 0.445 

1.148 

-0.670 

-0.309 


1.06l 

-.109 

-.190 

-m 

.045 

.046 

.145 

- 1.049 

.198 

-1.072 

-514 

-.843 

.301 

-1.277 

-583 

0.2b 2 

1-529 

-.283 

-.284 

.215 

.068 

.068 

.216 

- 1.504 

*301 

-1.469 

-.762 

-I.198 

•457 

-1.886 

-.881 


1.971 

-377 

-379 

.286 

.091 

.091 

.288 

-I.928 

.406 

- 1.833 

-.808 

- 1.524 

.6 17 

- 2.453 

- 1.152 


2.392 

-.473 

-.473 

•359 

.114 

.11 k 

.359 

- 2.329 

.513 

- 2.239 

-1.001 

-1.827 

.783 

-3.006 

-1.417 


2.802 

-570 

-.568 

.434 

• 137 

.136 

.431 

- 2.714 

.623 

- 2.383 

-1.139 

-1 

- 2.114 

•954 

- 3.554 

- 1.684 


0.689 

-0.226 

-0.210 

0.177 

0.056 


0.154 

-0.66^ 

0.140 

- 0.645 

-0.262 

-0.552 

0.185 

-O.85O 

-O.403 


1.244 

-.481 

-.466 

.370 

.H 7 


■ 349 

- 1.179 

.276 

-I.O77 

-.452 

-.940 

.422 

-1.600 

-.762 

0 . 4 b 2 

1-734 

-.717 

-.669 

.561 

.177 

.153 

.492 

- 1.613 

.424 

-1.562 

-501 

- 1.254 

.674 

- 2.311 

-1.102 

“SE - 

2.252 

-.961 

-.832 

.768 

.243 

.393 

.609 

-2.038 

.530 

-I.5H 

-.417 

-1.588 

.927 

- 3.115 

- 1.509 


2.858 

-1.212 

-.988 

■999 

.316 

.813 

.673, 

- 2.547 

.810 

- 1.530 

-.213 

- 1.983 

1.173 

- 4.090 

-I.966 


3-584 

- 1.470 

-1.186 

1.217 

.385 

.253 


- 3.140 

.961 

- 1.581 

•075 

-2.560 

1.635 

- 5.214 

- 2.546 
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TABLE 8 

CONVERGENCE OF ADDITIONAL CENTER DEFLECTION FOR DOUBLf CURVED 
SHEET UNDER NORMAL EDGE LOADS WITH 1/R Q MUCH 
GREATER THAN MERIDIAN CURVATURE 

(a) w 0 = h; Rq = 0 . 4 t> 2 /jdi 


— .2 

One term 

Two terms 

Four terms 
(a) 

2 2 
Ert h 

1.194 

-0.177 

-0.239 

-0.238 

2.077 

-.296 

-.400 

-.407 

3.074 

-.420 

-552 

-.562 


hj Rq = 0.2b 2 /jth 


Eit 2 h 2 


One term 


Two terms 


Four terms 


1.053 

1.932 

2.709 


-O.O83 

-.150 

-.208 


- 0.111 

-.205 

-.288 


- 0.111 

-.205 

-.288 


interpolated from four-term numerical calculations 
p^t ) 2 

for same values of — — — . 
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Figure 6 . - Additional deflection far sheets with i/Rq much greater 

than meridian curvature. 
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